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Abstract 

We characterize the avoided crossings in a two-parameter, time-periodic 
system which has been the basis for a wide variety of experiments. By 
studying these avoided crossings in the near-integrable regime, we are 
able to determine scaling laws for the dependence of their characteristic 
features on the non-integrability parameter. As an application of these 
results, the influence of avoided crossings on dynamical tunneling is de- 
scribed and applied to the recent realization of multiple-state tunneling 
in an experimental system. 



1 Introduction 

Avoided crossings of eigenvalue curves are generic features of quantum systems 
with non-integrable classical counterparts 1]. Their appearance allows for a 
wide variety of interesting, purely quantum mechanical, phenomena includin 
chaos-assisted tunneling |2(, the adiabatic exchange of eigenstate character 
and generally provides the mechanism by which underlying classical chaos affects 
the dynamics of a quantum system Their existence is also responsible for 
perhaps the most well-known result in the field of quantum chaos, the non- 
Poisson statistical distribution of level spacings in "chaotic" quantum systems 
0. 

In systems with two parameters, an avoided crossing along any curve in pa- 
rameter space can be associated to a "diabolical point" at which two eigenvalue 
surfaces become degenerate . The conical shape of the two eigenvalue surfaces 
in the vicinity of such a diabolical point ensures the characteristic hyperbolic 
behavior of two eigenvalues along any curve in parameter space passing near, 
but not through, the diabolical point. In the particular case of a near-integrable 
system, one parameter may be fixed to be zero leaving the system integrable for 
all values of the other parameter. Eigencurves will freely cross under variation 
of the latter parameter, thus creating diabolical points of the associated eigen- 
value surfaces when viewed in the full two-parameter space. As we show here, 
this type of diabolical point is important because perturbation theory can be 
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applied to characterize the conical shape and therefore characterize the avoided 
crossings of the near-integrable system. 

In this paper we study the particular two-parameter, near-integrable system 
of a harmonically driven pendulum: 

H(k, A) = p 2 + k cos 6 + A [cos (8 + tot) + cos (6 - vi)] . (1) 

This "one-and-a-half" degree-of-freedom system is one of the simplest types of 
classical systems to exhibit chaos. It is of significant experimental interest in 
quantum mechanics since it has been implemented in a number of studies 0, S 
IE Hfl LLI through the use of cold atom optics, particularly in investigations of 
multiple-state dynamical tunneling 0] • Theoretically, it provides a convenient 
framework for studying the avoided crossings of near-integrable systems since 
for A — > the system is the integrable pendulum Hamiltonian. 

In the following, we study the properties of avoided crossings for the driven 
pendulum with the use of Floquet theory. An avoided crossing of two Floquct 
eigenvalue curves (for A ^ 0) can be associated to a level crossing of the in- 
tegrable pendulum (A = 0) system and is characterized by the dependence of 
its closest approach on the non-integrability parameter A. For small values of 
A, we find that the spacing exhibits a power law dependence with an integer 
power. A modified degenerate perturbation theory is then applied to verify this 
dependence and associate it to the direct or indirect coupling of the associated 
integrable eigenstates. We then use the perturbation results to elucidate a mul- 
tiple state dynamical tunneling process in the vicinity of an avoided crossing 
and apply the results to the particular achievement of this tunneling in an atom 
optics experiment. We finally show the association of this avoided crossing to a 
nearby diabolical point. 

In Section [2 we present the model Hamiltonian under consideration in the 
paper, including a description of the system's classical dynamics. Section |21 
presents the quantum dynamics of the model system, with a brief review of Flo- 
quet analysis. Avoided crossings of the model system are investigated in detail in 
Section first numerically, then with the perturbation theory results presented 
in Appendix [XJ We review the implications of avoided crossings on dynamical 
tunneling in Section [S] and demonstrate the origin of those avoided crossings in 
a particular experimental system. Section contains some concluding remarks. 



2 The Model Hamiltonian 

The Hamiltonian we consider consists of a particle moving in the presence of a 
harmonically-modulated, spatially-periodic potential. It can be written in the 
form 

v' 2 

H'(p', x, t') = — + Vi cos(kx) + V 2 cos(kx) cos(wY), (2) 
2m 

where p' is the momentum and x the position of a particle of mass m, t' is time, 
V\ is the amplitude of the spatially periodic potential, V2 is the amplitude of the 
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modulation potential and oj' is the frequency of the modulation potential. The 
experimental implementation of quantum systems of with this type of Hamil- 
tonian was first proposed by Graham, Schlautmann, and Zoller in 1992 Q and 
then achieved by Raizen et, al. @, 0, , and Hensinger et. al. . 

It is useful to change to dimensionless units. We define: p = fr, 9 = kx, 

t = H f , co = k = A = ^, and H = g , where £ = Then, 

the Hamiltonian in Eq. Stakes the form 

H(p,6,t) = H [p, 6) + A [cos(0 - wt) + cos(6> + tot)] , (3) 

where 

H (p,9) =p 2 + K cos(6), (4) 

is the Hamiltonian of a pendulum and we have written the modulation term 
explicitly as two travelling waves. Note that momentum is measured in units of 
hk. 




Figure 1: Strobe plots of the system in Eq. with u> = 24 and parameters: 
(a) k = 7.8 and A = 0; and (b) k = 7.8 and A = 1.0. 



The classical phase space of a time-periodic one-and-a-half degree-of-freedom 
system such as H(p, 9, t) can be visualized by a strobe plot of the trajectories at 
times t = v— (y G Z), A strobe plot of phase space trajectories for the system 
governed by Hamiltonian H (p, 9, t) is shown in Figure^ a with parameters re = 
7.8 and A = 0. (For the case of A = 0, the system is independent of time and 
could be visualized with an ordinary parametric plot of phase space, however 
we plot the strobed phase space for convenience of comparison to the perturbed 
system). Because this system is integrable, all orbits lie on tori (in either the 
regions of the pendulum's libration or rotation) and the phase space is absent 
of chaos. Figure shows a strobe plot of the phase space with parameters 
re = 7.8, A = 0.5, and lu = 24. The travelling waves in the modulation term have 
phase velocities v — ±uj and are seen as primary resonance structures at p = ± j 
where 9 = v. Although much of the orbit structure of the integrable system is 
preserved, the tori with rational winding numbers have been destroyed, giving 
rise to a self-similar set of daughter resonance structures (see, for example, the 
two-island chains at p — ±j). Regions of chaos surround these resonances, 
most visibly near the separatrix of the pendulum resonance. 

3 Quantum Dynamics 

The dimensionless Schrodinger equation for the system in Eq. J3J is i^\>p{t)) — 
H(t)\ip(t)), where 

H(t) = p 2 + re cos 9 + 2A cos Uty cos(wi) , (5) 

We will consider the configuration space, 9 £ [0, 2n), to be periodic such that 
(9 + 2ir\ip(t)) = (9\ip(t)) and the momentum operator has integer eigenvalues: 
p\p) = n\p) (n £ Z). In the experimental systems, this is approximately 
achieved naturally because momentum transfer occurs in discrete units of hk 

When A = 0, the Hamiltonian reduces to that of the quantum pendulum, 
Hq(k) = p 2 + K cos 9 . The eigenstates of Hq are Mathieu functions 01 which 
we will henceforth denote as |"-(k)) so that _ffo(«)|n(«;)) = E n {K)\n(n)) where 
n = 0, ±1, ±2, . . .. (We will suppress the K-dependence of these eigenstates until 
their specification is necessary). States \n) with positive integer labels are those 
with even parity, states with negative integer label are those with odd parity. 
Note that as k — > 0, E n — > n 2 . If re ^ 0, but \n\ is large (i.e. the corresponding 
classical pendulum energy is much larger than that of the separatrix), we will 
again have E n w n 2 . 

3.1 Floquet Theory 

The Hamiltonian in Eq. © is time-periodic and therefore Floquet 's theorem 
guarantees that solutions of the Schrodinger equation can be written in the form 
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|1M*)> = e^Uft)) with |^a(* + T)> = K(t)), (6) 
where we have defined T = 2j and where \<j> a (t)) and f2 Q are called the Flo- 
quet eigenstate and eigenvalue, respectively. Substituting this solution into the 
Schrodinger equation yields the eigenvalue equation 

H F {t)\<t> a {t)) = \H(t) - i^j \<j> a (t)) = V a \Mt)) , (7) 

where Hp(t) is called the Floquet Hamiltonian. 

The Floquet Hamiltonian is a Hermitian operator in a composite Hilbert 
space Q<E)T [lilllfij . where 6 is the space of all square-integrable functions f{&) 
on the configuration space and T is the space of all time-periodic functions a(t) 
with period T and finite \ a {t)\ 2 dt. The inner product of two vectors \<p a ) 

and 1 4>b) in this space is then defined by 

= ^ (4> a \t){t\cl> b )dt = i / (Mt)\Mt))dt, (8) 

1 J-T/2 1 J-T/2 

where (</> a (i)|0&(i)) = {<f> a (t)\9){0\<f>i,{t))d$ is the usual inner product in 9. 
We select a complete orthonormal basis in this composite space 

(t\n,q) = \n)e^ t (n 7 qeZ), (9) 

where {\n}} are the eigenstates of the pendulum Hamiltonian Hq. These basis 
vectors satisfy ((n,q\n' ,q')) = S nin i5 q>q t. 

The Floquet Hamiltonian Hp is Hermitian, so the Floquet eigenvalues £l a 
are real and two Floquet eigenstates \(j> a ) and \(f>p) belonging to different eigen- 
values are orthogonal. Additionally, the Floquet Hamiltonian commutes with 
the parity operator defined by its action on the momentum eigenket n|p) = | —p) . 
Therefore the two operators can be diagonalized simultaneously and all Floquet 
eigenstates have definite parity: TL\<j> a ) = ±l\(f> a ). Floquet states with parity 
eigenvalue +1 will be called even, those with eigenvalue —1 odd. 

Given one Floquet eigenstate \(/> a (t)) with Floquet eigenvalue Q a , there will 
be another Floquet eigenstate \<t>' a (t)) such that \</)' a (t)) = e lqult \(t) a {t)) (q E 
Z), with eigenvalue Q! a = Vl a + qui. These two Floquet eigenstates, however, 
represent the same physical state, i.e. 

e- 4f °KW> = e-^UW) ■ (10) 

Therefore we may limit consideration to the fundamental zone —lu/2 < < ui/2 
in which each physical eigenstate of the time-dependent Schrodinger equation 
is represented by the corresponding Floquet eigenstate with eigenvalue within 
that range. 

Consider the unperturbed system H F = Hq — ijj which we will call the 
Floquet pendulum. The eigenstates of this system are precisely the basis states 
|n, q) with eigenvalues Cl nq = E n + quj. Figure shows the lowest nine energies 
of the even-parity eigenstates of the quantum pendulum and the corresponding 
Floquet eigenvalues in the fundamental zone —lu/2 < < lu/2. 
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Figure 2: (a) Energy curves of the nine lowest-energy, even-parity eigenstates of 
the quantum pendulum, H — p 2 + ncos9. (b) The nine corresponding "Floquet 
pendulum" eigenvalues in the fundamental zone — ^ < < +^ with u> = 24. 
The labels (n, q) on each Floquet eigenvalue segment identify the corresponding 
Floquet eigenstate |n, q). The dashed lines in (a) indicate the Floquet eigenvalue 
zone boundaries. 
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3.2 Another method for determining Floquet states 

An arbitrary dynamical state of the system can be expanded, with the use of 
equation ©, in the basis of Floquet eigenstates, 

a 

where the "prime" indicates that the sum is restricted to those Floquet states 
with Q a in the fundamental zone. The expansion coefficients are independent 
of time and can be written A a — (cj) a (0)\ip(Q)}. Using the time-periodicity of 
the Floquet eigenstates, we can then write 

hA(T)) =^e^° T |0 Q (O))(^(O)|^O)) =U(TM0)), (12) 

a 

showing that the time-evolution operator over a single period T 

U(r) = £V in « T |0 Q (O))(^(O)| (13) 

a 

is diagonalized by the Floquet eigenstates at time t = 0. We can therefore 
determine these time-strobed Floquet states by constructing the matrix U TO m' = 
(m|U(T)|m') in some convenient basis {|ttt,)} in O, truncating this matrix at 
some appropriate level m — M where it becomes approximately diagonal (i.e. 
Umm >> Umto for m ^ M), and then performing a numerical diagonalization 
to obtain the \<p a {0)) and Sl a (mod u). The m th column of U is obtained by 
evolving the basis vector \m) over one period T via numerical integration of the 
Schrodinger equation. 

In subsequent sections, we will compare the phase space distributions of the 
time-strobed Floquet eigenstates |</> Q (0)) to the classical system. We can do this 
by introducing the Husimi distribution p(0o,po) (l6lll7| of a quantum state \<j>) 
on the classical phase space (0o,Po) 

p(9 , Po ) = -L|(0 O7PO |0)|2 ; (14) 

where the coherent state \0q,Po) is defined as an eigenstate of the annihilation 

operator a — (o/o- + iaf^j with position and momentum expectation values 

of 8q and po respectively. The free parameter a is set according to the physical 
system considered (see below). The representation of such a coherent state in 
the discrete momentum basis {\p}} is given by 



(p\0 o ,Po) = Aexp 



(7^ 

y(p- Po) 2 - i0o{p-Po) 



(15) 



where A is a normalization factor guaranteeing (0o,Po\0OtPo) = 1- The action of 
the annihilation operator on the coherent state can be used to show that (p) = 
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(e ,Po\p\e ,Po) = Po, (0) = , A8 = a/V2, and Ap = (a^)" 1 . Thus, the 
coherent state is a minimum- uncertainty wavepacket, where the free parameter 
determines the ratio of its uncertainty in position and momentum, i.e. a 1 — 
AO/ Ap. Reference 01 presents an in-depth discussion on the selection of the 
parameter a. In all Husimi plots shown in subsequent sections, we set a — 
1.18 k -1 / 4 , a choice which provides the best association between the quantum 
pendulum eigenfunctions and the corresponding classical orbits). As we will see, 
the Husimi distributions of the Floquet states lie directly on the orbit structures 
of the classical phase space. 

4 Avoided crossings 

The Floquet pendulum is integrable and its eigenvalues U, n q {ti), shown in Fig- 
ure|21b, cross under the variation of n. For any nonzero A, however, the system 
represented by the Hamiltonian in Eq. (jSJ is non-integrable and the approach 
of any two (same parity) Floquet eigenvalues under variation of k results in 
an avoided crossing. This well-known result, the no-crossing theorem, was first 
proven by von Neumann and Wigner for eigenvalues of generic Hermitian matri- 
ces Q . They also showed that adiabatic passage of two quantum states through 
an avoided crossing leads to an exchange of character. (In a two-parameter 
system, this exchange can be related to the partial circuit of a diabolical point 
6] , while in single parameter systems it can be related to exceptional points in 
the complex parameter plane |l8j). Avoided crossings of Floquet eigenvalues in 
the fundamental zone, which generally involve states localized in well-separated 
regions of the phase space, will therefore allow a wide variety of interesting 
quantum dynamical phenonomena, including adiabatic transitions and tunnel- 
ing. 

In this section, we will consider the near-integrable regime (0 < A << k), in 
which a clear association can be made between the Floquet eigenstates of the 
perturbed system (A ^ 0) and those of the Floquet pendulum (A = 0). In this 
regime, the Floquet eigenvalues will follow nearly the same dependence on k 
as the unperturbed eigenvalues seen in Figure [21b, except in the vicinity of an 
avoided crossing. For k values sufficiently far from these avoided crossings, we 
can make a unique, though necessarily local, association \4> a ) «-> \n a ,q a ) of the 
Floquet eigenstate \<j) a ) to the Floquet pendulum state with maximum overlap 
|(n, q\(j> a )\- As A — + 0, this association will become an equality. We will see that 
the fundamental characteristics of an avoided crossing between states \4> a ) and 
will be determined by the difference 

Aq al 3 = \q a - qp\ . (16) 

In the subsections below, we first present a numerical analysis of some represen- 
tative avoided crossings in the system with u> = 24, and then use perturbation 
theory to show that the results are quite general. 
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Figure 3: A Aq a p — 1 avoided crossing of the system Hp with parameters 
lo = 24 and A = 5 x 10~ 2 . The Husimi distributions of the corresponding 
Floquct eigenstates are shown at K = (7.455, 7.460, 7.465) (the horizontal axis 
is 9o & [0, 27r), vertical is po € [—15, 15]). The dotted lines are the eigencurves 
of the unperturbed Floquet pendulum. 



4.1 Numerical Results 

Three avoided crossings of Floquet eigenvalues in the fundamental zone with A = 
5 x 10 -2 are shown in Figures|310J and[5]with Husimi plots of the corresponding 
Floquet eigenstates overplotted on the figures. The dotted lines shown are the 
eigenvalues of the unperturbed Floquet pendulum. These plots were created 
by numerically calculating the time-strobed Floquet states at a sequence of n 
values. With each step forward in k, the new states were associated to those of 
the previous step by calculating the maximum overlap and verifying continuity 
of the eigenvalues. In the case that two Floquet eigenvalues crossed between k- 
steps, the size of the step was reduced and the process repeated until no crossing 
occurred. 

Each of these three avoided crossings involves one Floquet eigenstate local- 
ized within the pendulum resonance at p = and another localized outside of 
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Figure 4: A Aq a p — 2 avoided crossing of the system Hp with the same A and 
u values as in Figure [3] The Husimi distributions shown are the two Floquet 
eigenstates at k = (7.82146,7.82150,7.82154) (axes are the same as FigureEJ). 
The dotted lines are the eigencurves of the unperturbed Floquet pendulum. 
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Figure 5: A Ag Q/ 3 = 3 avoided crossing of the system Hp with the same A and 
uj values as in Figure [3] The Husimi distributions shown are the two Floquet 
eigenstates at k = (9.973221751, 9.9732217534, 9.973221755) (axes are the same 
as Figure EJ. The crossing of the corresponding eigencurves of the unperturbed 
Floquet pendulum falls outside of the plotted region at k sa 9.973242. 
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the pendulum resonance. The avoided crossing in Figure |3] involves the states 
(n a> q a ) — (1)0) and (np,qp) = (5,-1); Figure 0] involves the states (n a ,q a ) — 
(1,0) and (np 7 qp) — (7,-2); and Figure involves states (n a ,q a ) = (0,0) 
and (np,qp) — (8,-3). The associated crossings can be found in Figure Elb. 
These particular avoided crossings were chosen as representative examples with 
Aq a p = 1, 2, and 3, respectively. 

Some general characteristics of these avoided crossings deserve attention. 
First, the "exchange of character" between the two states is evident in the 
evolution of the Husimi distributions with k. The associations \4> a ) \n a ,q a ) 
and \4>p) «-> \rip,qp) well before the avoided crossing become \<f> a ) <-> \np,qp) 
and \4>p) <-> \n a ,q a ) well after. For k values at the avoided crossing, the two 
Floquet states are superpositions of the asymptotic states. Second, there is 
quite a disparity of scale among the three avoided crossings. In particular, the 
minimum eigenvalue spacing A a p, defined by 

A a p =min(|0 Q -£lp\) , (17) 

is relatively large for the first and relatively small for the third. Finally, the 
position K ac of the minimum spacing (which we will henceforth call the "position 
of the avoided crossing"), is not necessarily equal to the position kq of the 
unperturbed crossing. Indeed, it seems that for Aq a p ^ 1, the avoided crossing 
is significantly offset in both k and tt. 

To make these last two observations more quantitative, we have computed 
dependence of A a p and An ac = \kq — n ac \ on the parameter A. The results for 
the three example avoided crossings are shown in Figures El and We see that, 
for small values of A, the dependences are all well approximated by power laws 
with integer exponents. The minimum spacing of the avoided crossings is given 
by 

A a p = A X A ^ , (18) 

where the coefficients are A w {5.44 x 10~ 3 , 1.02 x 10~ 3 , 7.35 x 10~ 7 } for avoided 
crossings I, II, and III, respectively. The K-offset of the avoided crossings are 
given by 

A Kac = B A 2 , (19) 

where B « {5.2 x 10 -3 , —8.3 x 10~ 3 } for avoided crossings II and III (Ak„ c = 
for avoided crossing I). 

4.2 Perturbation theory results 

We now use perturbation theory to determine the behavior of the Floquet eigen- 
values and eigenstates in the neighborhood of avoided crossings. We will obtain 
approximate solutions (\4> a ), Q a ) to the Floquet eigenvalue equation 

Hp(K,X)\4> a (K,X)) = |0a(«,A)) = n a (« l A)|0„(K J A)) 1 (20) 

where ffp(fc) is the Floquet pendulum Hamiltonian, V — 2cos(? cos(wi), and 
A is considered a small expansion parameter. Our unperturbed system is the 
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Figure 6: The minimum spacing A a p as a function of A of the three avoided 
crossings shown in Figures |3 0] and[S] The functions A Q ^ = 5.44 x 10 _3 A, 
A Q/3 = 1.02 x 1CT 3 A 2 , A a0 = 7.35 x 10" 7 A 3 are overplotted. 
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Figure 7: The /c-offset AK ac as a function of A for the avoided crossings shown 
in Figures E| and The functions |AK ac | = 5.2 x 1CT 3 A 2 and |AK ac | = 8.3 x 
10~ 3 A 2 are overplotted. 
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two-fold degenerate system Hp(Ko), where n is the parameter value at which 
the eigenvalue curves of two Floquet pendulum states \a (k)) = \n a ,q a ) and 
\/3°(k)} = \np, qp) cross. We have seen in Section ll~Tl that, for A 7^ 0, the closest 
approach of the eigenvalues O a and £lp involved in an avoided crossing may not 
occur at k = kq, so an offset must be allowed for. We therefore introduce into 
l(2T?|l an arbitrary function n(\) = n + Ak(A) and expand Ak(A) as a power 
series in A. The particular value K ac (A) at which the eigenvalues make their 
closest approach can then be determined by solving the extremal condition for 
AQ a p = |0 o (/c, A) - top(K, A)|, 

dAn a p 



da 



0, (21) 



at each order to fix the expansion coefficients of Ak(A). In this manner, we find 
the perturbed eigenstates and eigenvalues at k = n ac . 

The details of the perturbation analysis are given in Appendix 1X1 The re- 
sults may be summarized as follows. The breaking of the degeneracy between 
states I a ) and \/3°) occurs at the lowest order X N for which the coupling be- 
tween these states, v^l (defined below), is nonzero. At this order, the two 
Floquet eigenvalues in the region of the avoided crossing are determined by the 
eigenproblem of a 2 x 2 matrix in the basis of the unperturbed states \a°(no)} 
and |/9°(ko)>: 



A K W6E a + v { a V ig> 




(22) 



V vfj An^SEp 

The coefficients and determine the zeroth-order near-degenerate eigen- 
states in the region of the avoided crossing; Ak' ' and Of 5 are the coefficients 
of X N in the expansions of the arbitrary K-offset and the near-degenerate Flo- 
quet eigenvalues fi±, respectively; 5Ei = (nj(«o)| cos0|rij(/co)) are the slopes of 

(N) 

perturbation operator: 

Vim = ((nuQl\V\n m ,q m )) = (m(Ko)\ cos 6\n m (K )) {S qi ,q m +i + <W m -i) ■ (23) 
For the first three orders in A, these couplings are 



the unperturbed eigencurves; and depends on the matrix elements of the 



(i) 

ij - "ij 



= Vn (24) 



v (2) _ \ " ViyVyj 



13 ^ 0(0) _ o (0) 

y (3) _ \- ) \ - VyyVytrVaj 



(25) 




a ^p } (oW-fif) (0(0) 



, x Vi-yVyi [V aa (SE a - 6E 7 ) + Vpp (8E 1 - 6Ep)} l 



o(°) - a, 



(0) 
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™ac 


(n a ,q a ) 


(np,qp) 


Aq a p 


A 

n num 


A p t 


^num 


B p t 


7.46 


(1,0) 


(5,-1) 


1 


5.44 x 10" 3 


5.5 x 10" 3 






7.83 


(1,0) 


(7,-2) 


2 


1.02 x 10~ 3 


1.0 x 10~ 3 


5.2 x 10~ 3 


5.2 x 10~ 3 


9.97 


(0,0) 


(8,-3) 


3 


7.35 x 10~ 7 


7.3 x 10" 7 


-8.3 x 10" 3 


-8.3 x 10" 3 



Table 1: Comparison of numerical results to those of perturbation theory. The 
quantities A num and B num are obtained from numerical simulation and A pi and 
Bpt are obtained from perturbation theory. 

where i,j € {a, (3}, the zeroth-order eigenvalues Q ; are taken at kq, and 
we write rj (0) = fl^ . Using Eqs. (J22J) and QU we find that the N th -order 
corrections to the eigenvalues, at the position of the avoided crossing, are given 
by 



( K nc) — 2 



± 



} aa T Vpp 



+ Ak^ } (5# a - SE P ) 



+4i«iyi 2 



5£ a 



(TV) 

^ 1 ^ 1 



(27) 



At orders < M < N, the two eigenvalues tt± are degenerate at an offset from 
(«o,f^°)) specified by the coefficients 



(M) _ CM) 
Upp Vaa 

SE a - SEp 



(28) 



and 



j pp 



SE a 



Vaa Otip 



5E a - SE, 



'0 



(29) 



The origin of the numerical results presented in section 14.11 is now clear. 
The matrix elements of the perturbation Vij are non-zero only when Aq^ = 1. 
Therefore, for an avoided crossing between states \(j> a } and \<j>p), we must have 
Aq a0 = N. Using Eq. Q2J and the fact that f)( M<Ar ) is the same for \<f> a ) and 
\<f>p), we find that, to lowest order in A, the minimum spacing is given by 



la/3 



„(A 9q/3 ) 
J af3 



(30) 



Its interesting to note that the K-offset of an avoided crossing in this system is 
dependent on A 2 , because Val and Vgp are non-zero, even when Aq a p > 2. 

Table ^ shows a quantitative comparison of the numerical results presented 
in the previous section (in terms of the coefficients A and B of Eqs. I|18l) and 
l|19[)) to those obtained by the perturbation analysis. For all three examples 
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of avoided crossings, we see excellent agreement. We have also verified the 
predictions of perturbation theory for a number of other avoided crossings in 
this system (as well as those in systems with different values of oj), finding 
similar agreement. 

A number of other characteristics of the avoided crossings can be determined 



from our perturbation analysis. Substituting fi± (/«ac) and Akuc' into Eq. 
I|22|). we find that at the position of the avoided crossing, the two perturbed 
Floquet eigenstates (to lowest order in A) become an equal superposition of the 
two associated Floquet pendulum states, i.e. 



\Ct(A Kac 



\C±(Ak 



(31) 



We may also determine the relative magnitudes of these coefficients at some k 
value near the avoided crossing. If, instead of calculating Ak' w ' by the extremal 
condition, we instead determine the K-offset where the eigenvalue separation is 
a times the minimum value, we find 



AN) 



± 



yj3 SE aa - SEp/s \SE aa - 6Epp\ 
A simple calculation then shows that the coefficients obey 



\C±\ 

\cf\ 



1 



iA a 



(32) 



(33) 



where we have assumed that (6E a — SEp) > 0. As a — ► (0, oo), we see that, for 
example, j-^j — > (1,0), as expected. This verifies the qualitative behavior seen 
in the Husimi distributions plotted in Figures 13151 



5 Implications for Dynamical Tunnelling 

Under the evolution of the system given in Eq. JSJ, an initial quantum state 
|'0+(O)) which is localized in classical phase space (in the sense of its Husimi 
distribution) in a region of positive momentum at p ss po may undergo time- 
periodic dynamical tunneling, across the central resonance (and all intervening 
KAM tori) to the opposite momentum region at p —po- The mechanism 
for this behavior is the existence of a near-degenerate and opposite parity pair 
of Floquet eigenstates which each have localization near po and — pq. If the 
Floquet eigenvalues of these two states are far from any avoided crossings, the 
tunneling dynamics are well described by a two-state process exactly analogous 
to the tunneling through a potential barrier in the time-independent double well 
system |20j . In the vicinity of an avoided crossing, however, the dynamics are 
influenced by a third state with partial localization in the regions of ±po and the 
time-evolution takes on a more complicated beating behavior. In this section 
we analyze this tunneling behavior in the perturbative regime (A small) and 
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Figure 8: The avoided crossing from Figure 01 (a), now with the relevant odd- 
parity state included (dotted line). The vertical lines are the n values at which 
the time-evolution in Figure [!|] was obtained. 

then apply the results to tunneling oscillations observed in a recent experiment. 
Although the experimental system cannot be considered to be in a perturbative 
regime, we identify the diabolical point associated to the relevant avoided cross- 
ing and show that an approximate result can be obtained numerically which 
characterizes this avoided crossing quite well. 

5.1 Tunneling behavior in the perturbative limit 

To analyze the tunneling induced by the Hamiltonian in Eq. JSJ), we consider 
again the one-period time-evolution operator U(T) defined in Sect ion and its 
eigenvectors, the time-strobed Floquet eigenstates \4> a {T)) — |0 a (O)) (we will 
drop the explicit reference to "time-strobed" in this section). As a particular 
example, consider the avoided crossing shown in Figure |H1 (the same as that 
shown in Figure but with the odd-parity state now included, shown as a 
dotted line). At k = 7.456, a value far from the avoided crossing, the opposite- 
parity pair of states \4> a ) and |</> 7 ) are near-degenerate (|fi 7 — fi Q | w 4.5 x 10~ 6 ) 
and have localization at p ±5. The phase space localization of state \4>p) is 
completely within the central resonance and does not overlap significantly with 
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this pair. We can construct an initial state, localized at either p 5 or —5, as 
an equal superposition of the two near-degenerate Floquet states, i.e. ^±(0)) = 

(\4> a ) ± l^))- When either is acted on by U(nT), the evolution is periodic, 
oscillating between p ss ±5 with a tunneling frequency uotun = — fi 7 |. The 
corresponding number of modulation periods for complete oscillation is then 
n tU n = u/ujtun ~ 5.3 x 10 6 . 

In the region of the avoided crossing, the odd-parity state |<^ 7 ) will be ap- 
proximately unchanged, while states \4> a ) and \4>p) become, to lowest order in A, 
superpositions of their unperturbed, Floquet pendulum, counterparts. There- 
fore, in the neighborhood of the avoided crossing, both \(j> a ) and \<fip) will have 
significant support in the same region of phase space as |0 7 ). At the exact posi- 
tion K ac of the avoided crossing, where \4> a ) and \4>ii) are an equal superposition 
of the unperturbed states, the initial conditions localized at either p w 5 or — 5 
can then be written 

\M0)) = \\<t>*) + l\<t>p)±A=\<h), (34) 

where all three eigenstates are evaluated at k = n ac . Applying the time- 
evolution matrix to |^+(0)) we find, after n applications, 

|V+(nT)> = U(nT)|V+(0)) (35) 

' e -mTAn Q/3 e -inTAfi 7/3 
— \<l> a ) + gl^) + 10 

where Af2y = fij — Slj, here. If we make the assumption that |AQg T | = 
|Af2 a7 | = A at n = n ac , we see that ~ IV'-(O)). This three-state 

process therefore generates a new, larger, tunneling frequency uj tU n = A. The 
time-evolution of \ip+(0)} at k k n ac is shown in Figure 0a. Notice that the 
theoretical tunneling period n tun — ^ ~ 1.8 x 10 5 is modulated by a beat period 
resulting from the small difference A0^ 7 — Af2 7Q . As the position of the avoided 
crossing is chosen more precisely, and i7 7 — > Qc "~^ n ' 1 ; the period of beating goes 
to infinity. 

Between these two extremes of regular two-state and three-state tunneling, 
at other values of the parameter k along the avoided crossing, the tunneling 
takes on a beating behavior due to the two eigenvalue differences between the 
odd-parity state and each even-parity state. An example of this (the evolution 
of |^+(0)) under U(T) at K = 7.4605) is shown in Figure 0b. For this and 
all values of the parameter k, the evolution of the momentum expectation of 
\ijj+(0)) is well fit by the simple function 

(p) (nT) = A a cos (An ja nT) + Ap cos (Aft^nT) , (36) 

where A a and Ap can be related to the overlap of \ip+(0)) with \<fi a ) and \4>/3), 
respectively. 

The variation of two-state tunneling frequencies in the vicinity of avoided 
crossings has been remarked on by many authors 0, 121I E2I I23I E3. 125L [2(| in 
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Figure 9: The evolution of state |^ + (0)) (created by the superposition of the 
near-degenerate states at k = 7.456) under the action of U(T) at n m n ac (a) 
and k = 7.4605 (b). 
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many different systems, and is often attributed to the influence of underlying 
classical chaos. Classical chaos in a non-perturbative regime will certainly in- 
troduce additional complications to the quantum dynamical tunneling process 
which we have not investigated here (notably the interaction between tunneling 
through dynamical barriers and free evolution in a region dominated by chaos 
|27|^: and avoided crossings will become larger, more numerous, and may over- 
lap, leading to interaction of more than three relevant states. However, as we 
have seen in this section, the basic mechanism for "tunneling enhancement" 
does not necessarily require global chaos but only the non-integrability which 
leads to avoided crossings. Indeed, for the non-integrability parameter used in 
this section (A = 5 x 1CP 2 ), the classical phase space has only small regions of 
chaos. 

5.2 Analysis of a Tunneling experiment 

As an application of the tunneling results of the previous section, we consider 
the experiment of Steck, Oskay and Raizen \l2\ . The Hamiltonian implemented 
in this atom-optics experiment depends on a single parameter, a, and can be 
considered a special case of that in Eq. (J5J by setting k sw a/2. 17 (up to a 
sign difference which can be removed by a 7r-translation of the angle variable), 
A = k/2, and u> « 6.04. It should be noted that this system is not connected to 
the Floquet pendulum system since A is not an independent parameter. Instead, 
in the limit A, k — > 0, the free particle Hamiltonian is obtained. 

The primary result of the experiments detailed in |12j is the observation of 
dynamical tunneling, between two resonance islands in the classical phase space 
(located at p ss ±3), which exhibits oscillation frequencies dependent on the 
parameter value considered and independent of the modulation frequency. In 
particular, for a range of parameter values, the observed tunneling oscillations 
were dominated by two primary frequencies. The Floquet eigenvalue curves 
for the experimental system are shown in Figure I1UI and the experimentally 
observed tunneling frequencies are shown in Figure llll b. overplotted on the 
differences of eigenvalues between three particular Floquet eigenstates. These 
three states exhibit significant localization in the region of the classical resonance 
islands at some values of the parameter k |ll| . The experimental frequencies 
are well predicted by only two of these differences, namely Af2 Q7 and AQpy. 

The oscillations seen in the parameter region k £ [3,4.75] are due to the 
existence of an avoided crossing between the eigenvalues of the even-parity states 
labelled \4> a ) and \4>p). As can be seen in the overplotted Husimi distributions, 
these two even-parity states originate at small k values from two disconnected 
regions of the phase space, the first residing at p « 0, and the other at the 
positions of the classical resonance islands where the odd parity state |</> 7 ) also 
has its primary support. In passing through this avoided crossing, the odd-parity 
state retains its original character, while the two even-parity states become 
mutual superpositions of their small n character. It is clear that this avoided 
crossing is not of the ideal form considered in the perturbative regime. The 
superposing effects of this avoided crossing are extend well beyond the minimum 
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Figure 10: Floquet eigenvalue curves for the system considered in G3 (A = k/2, 
uj » 6.04). 
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Figure 1 1 : The boxed region of Figure EI] where an avoided crossing occurs 
between two even parity states (a). The Husimi distributions of these and a 
third, odd-parity, state are overplotted, with axes of 9q £ (0, 27r) andp € (—6, 6). 
The differences of these three eigenvalues are shown in (c), with experimental 
tunneling frequencies overplotted (circles, reprinted with permission from Steck, 
et. al. Fig. 1). 



eigenvalue separation at re » 3.0 (see FigureEJa), due to the natural eigenvalue 
curvature of a state with the low re character of \<j>p) (a simpler example of this 
is state (1,0) in Figure|21b). Just past this minimum spacing, the eigenvalue of 
state \4> a ) curves back downward toward that of \4>p), immediately beginning a 
second avoided crossing and thus preserving the composite nature of these two 
even-parity states through re « 6. As a further complication, in the parameter 
region re £ [5, 7], states \(f> a ) and \<frp) are joined by two other even-parity states 
in a complex and overlapping set of avoided crossings. 

Despite these complications, we can identify the dynamical tunneling in the 
parameter region re £ [2.5, 4.75] to be a three-state process involving those states 
shown in Figure ITT1 a. As predicted by the results of the previous section, the 
observed tunneling frequencies involve only differences in Floquet eigenvalue 
between the odd parity state and the two even parity states. A direct com- 
parison of the oscillations from ^3] with the numerically calculated evolution 
of an initially localized state |^>+(0)) under U(T) is shown in Figure IT51 for the 
values re = 3.68 and 4.50. Neglecting dissipation and a momentum offset of 
the experimental values (due to the fact that not all atoms contributing to the 
average are participating in the dynamics), there is good agreement in the sec- 
ond case. In the first case, and for all parameter values between re = 3 and 4, 
the experiment seems to pick up only one of the underlying frequencies (Ar2 7/ g, 
while the numerics predict nearly equal contributions from Afi 7/ g and Af2 Q7 . 
It was noted in Reference [i^ that the detection of fewer than the predicted 
number of frequency contributions to the tunneling behavior was also found in 
another experimental system jjjj- 

Finally, we would like to consider the origin of the avoided crossing involved 
in the dynamical tunneling observed in [l2T |. Figure H^l shows the avoided cross- 
ing of Figure ITT1 a lying on the eigenvalue surfaces tt a and fl^ in re — A space. 
One can see that these two surfaces meet at a diabolical point on the re = 
axis, where A ~ 2.8. Numerical analysis shows that, in the neighborhood of the 
diabolical point, the minimum spacing between the two eigenvalue surfaces is 
linearly dependent on re, with 



This result is in good agreement with a perturbation analysis similar to that 
of Appendix A, but with re as the small expansion parameter. The predicted 
minimum spacing from such an analysis yields 



where the matrix element must be numerically calculated as the coupling of 
the two degenerate Floquet eigenstates of the (re = 0.0, A = 2.8) system through 
V' = cos 9, the coefficient of re. Although the avoided crossing in the experiment 
(Figure ^2 a) appears at a parameter value outside the range of validity of 
perturbation theory, the observed minimum spacing between ASl a and Aflp 
agrees with Eq. 13131 to within a factor of two. 



A a p = 9.22 x 10~ 2 re. 



(37) 



A a/3 , pt = 2|^ /3 |re«9.2x KT 2 re, 



(38) 
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Figure 12: Numerical evolution (squares) of a positive momentum centered 
initial state under the time-evolution operator at k = 3.68 (a) and 4.50 (b), 
with experimental values overplotted (circles, reprinted with permission from 
Steck, et. al. [13 , Fig. 2). 
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Figure 13: Eigenvalue surfaces, over n — A space, of the even parity Floquct 
cigcnstates involved in the avoided crossing of Figure lllla. The eigenvalues 
on the particular curve investigated in reference [l2j are overplotted in bold. 
A dashed line traces the minimum separation of the two surfaces from the 
diabolical point at k — 0, A ~ 2.8 to the avoided crossing seen in the experiment. 
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6 Conclusions 



We have shown that the minimum spacing of avoided crossings in a near- 
integrable time-periodic system, and therefore the minimally separated "ridges" 
of the double-cone structures surrounding a diabolical point, exhibits a power 
law dependence on the non-integrability parameter, with integer power. A mod- 
ified degenerate perturbation theory has allowed us to relate the coefficient of 
this dependence to the direct or indirect coupling of the two related unper- 
turbed Floquet states through the perturbation operator, and the integer power 
to the number of "photon energies" to connecting their related energy eigenval- 
ues. Moreover, the perturbation analysis predicts a qualitatively identical struc- 
ture for all avoided crossings which allows us to characterize generically their 
affect on dynamical tunneling. This description was applied to a particular 
avoided crossing generating multiple-state tunneling oscillations in an experi- 
mental system, and its connection to a nearby diabolical point was revealed. 
It is hoped that these results will provide guidance for the development of new 
experimental setups which intend to use avoided crossings in the realization of 
multiple state tunneling processes and adiabatic transitions, or the preparation 
of Schrodinger's cat type superposed states. 
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A Perturbation Theory for Avoided Crossings 

The Floquet pendulum system H F (k) has a two-fold degeneracy at k = kq where 
two eigenstates |a°) = \n a (Ko), q a ) and |/3°) = \np{Ko), qp) have eigenvalues 
^a(Ko) = ^/3( K o)- We will use a modified degenerate perturbation theory to lift 
this degeneracy at the /t-position of the resulting avoided crossing when A > 0. 
In light of the discussion in Section 14.21 we write the Floquet Hamiltonian as 

H F (k) = H% + SH F Ak(A) + A V, (39) 

where we now write H F = H f (kq). The operator 8Hp is — cos 8 with the 
simplifying assumption that this operator is diagonal in the unperturbed basis. 
This is equivalent to the assumption that the unperturbed eigenvalue curves are 
linear in the small region of k under consideration. We expand Ak(A), setting 
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Arc(O) = 0, so that 



Are(A) = AkWa + Ak< 2) A 2 



(40) 



We expand the near-degenerate eigenstates and eigenvalues in powers of A about 
their (k — kq, A = 0) values: 



0) = C a \a°) + C \f3°) + AI0 1 ) + A 2 |0 2 ) + • • • 

n = fi (0) + Att« + \ 2 n^ + ... 



(41) 



where H F \a°) = Q^\a°), H F \(3°) = fl^\0°) and the zcroth-order eigenstates 
have been assumed to be a superposition of the two degenerate unperturbed 
states. As is usual in degenerate perturbation theory, the lowest-order near- 
degenerate eigenstates and the corrections to their eigenvalues will be the eigen- 
vectors and eigenvalues of a 2 x 2 matrix in the basis of the two degenerate unper- 
turbed states. At this order, we will make the distinct assignments \(j>) — * \4>±) 
(corresponding to the solutions of the quadratic characteristic equation). The 
final expressions for C*±, Cj, |<4), and n£ >0) will depend on k through the 
particular choice of the Area's. The Floquet eigenvalue equation Q for the 
near-degenerate state now takes the form 



(H° F -Q^) (C a \a Q ) + C P \P }) 



+ A 



(h° f - n^M) + (akW 



SHp + v - n 



(i) 



+ A 2 



(H F - ft (0) )|</> 2 ) + {Ak^SHf + V- ftW)^ 1 ) 



+ (A K ^6H F - rt 2 >) (C a \a°) + C P \P Q )) 



= . (42) 



A.l First-order results 

At first order in A we have the eigenvalue equation 

{H F - fiW)^ 1 ) + {An^5H F +V- n^){C Q \a°) + C^|/3 )) = 0. 
If we now act on this equation with (a 1 and then with (/?° | , we obtain 



AkW(5£; q + v aa - v aP 



0. 



(43) 



(44) 



where Vij are the matrix elements of the perturbation in the basis of the degen- 
erate unperturbed eigenstates, 



SEi = {{ni,qi\5H F \m,qi}) = (n,| cos0|n,) , 
and we have used the identity 



(45) 



l = El^°l^ E K<7-y><n 7 ,g 7 |= E MM- (46) 

7 ("t^t) 



n,q— — oo 
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At this point, we must consider two possible cases: V a p ^ and V a p — 0. In 
the first case, a nontrivial solution for the Cj requires first-order corrections to 
the unperturbed Floquet eigenvalues 



(i) _ 



1 



V aa + Vpp + An (1 \5E a + SEp) 

(V aa - Vpp) + A^\8E a - 8Ep)\ ' + A\V a p\ 2 
We can then write the spacing between the two Floquet eigenvalues 



±1 
2 



(47) 



(V aa - Vp P ) + A K ^(SE a - 5Ep)\ + 4\V a p\ 2 j ' +0(A 2 ) 

. (48) 

In order to determine the particular value of Ak' 1 ' at which the minimum 
spacing occurs, we solve the extremal equation 



(dAn\ _ 1 / dAfl 
\ 8Ak J > ~ A \8AkW 



(m^l), 



(49) 



A,AK (m ' 



(50) 



finding that 

A (i) = Ym - y »» 

ac SE a - SEp 

If we substitute this into Eq. 147(1 . we obtain first order corrections to the 
eigenvalues of 

,(i) _ VppSE a — V aa 5Ep 



5E a - SEp 

Therefore, we find the minimum splitting to be 

A a p = 2X\V a p\ + 0(X 2 ) 



± \V af i\ 



(51) 



(52) 



In the second case (V a p — 0), we find that the nearest approach of the two 
Floquet eigenvalues is in fact a crossing (to first order in A). In this case we 



define A«ic to be the offset of this crossing, i.e 

,2 



An = 0(\ 2 ) at A«W = V PP V * 



SE a - SEp ' 

The coefficients Ci remain undetermined until the degeneracy is broken. 



(53) 



A. 2 Second-order results (V a p = 0) 

For the case V a p — 0, the zeroth order states must be determined from the 
second order equation, which takes the form 

(h° f - n^)\4> 2 ) + [Ak w sh f + v- fiW)!^ 1 ) 

+ (A K ^SH F -n^) (C a \a°)+Cp\f3°)) = 0, (54) 
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Following the same procedure as in first order, we obtain 

AK^5E a +v aa - v aB \ ( C a 

where we have used the first-order result 



AkWSeJ+v^ - Q< 2 ) ) ( % ) ~ ° ' (55) 



(7V) = ^^y (7 (56) 



with n\ u> = ( n 7 , qj\Hp\n-y, <7 7 ), and where we have defined 

7^{a,,3} " ~ iZ 7 

Again we must consider two cases: w Q /3 ^ and v a p = 0. In the first case, our 

(2) 

procedure for determining An a( ! is identical to that of first-order and we obtain 



AK ^^ SE a - V SE p and A «0 = ^ 2 M + O(\ 3 ). (58) 
In the case that v a p — 0, we find that AO = 0(A 3 ) at an offset of 

A. 3 Third-order results (V a p = and v a p = 0) 

If the conditions V a p — and v a p = are satisfied, we can attempt to lift the 
degeneracy at order A 3 . We obtain the following results 



= 0, 
(60) 



where we have defined 



^feU (n(°)-O 7 0) ) (fi(°)-fi< 0) ) 

and 

«.<= X! - — — ; , 2 (ie{«,/3», (62) 

and where we have used an expression for (7°|</> 2 ) (analogous to Equation I56fl 
determined from the second-order equation. 

These equations are nearly of the same form as those at first and second 
order. By the same procedure we determine 

A"£ = * W ^Ijg~ ga and A«/ 3 = 2A 3 K /3 | + 0(A 4 ), (63) 

for the case that X> a p ^ and A a p — 0(A 4 ) when K> a p = 0. 
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